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ABSTRACT 


The Colorado Center for Astrodynamics has completed a study comparing ocean surface 
winds derived from GPS bistatic measurements with QuikSCAT wind fields. We have 
also compiled an extensive database of the bistatic GPS flight data collected by NASA 
Langley Research Center over the last several years. The GPS data are augmented with 
coincident data from QuikSCAT, buoys, TOPEX, and ERS. 

Our findings and recommendations are included in the following articles, which are 
attached. 

1. “Development in Using GPS for Oceanographic Remote Sensing: Retrieval of 
Ocean Surface Wind Speed and Wind Direction,” A. Komjathy, M. Armatys, D. 
Masters, P. Axelrad, V.U. Zavorotny, and S.J. Katzberg, Proceedings of the ION 
2001 National Technical Meeting , Long Beach, CA, January 22-24, 2001 . 

2. “GPS-based Remote Sensing of Ocean-Surface Wind Speed from Space,” 
Armatys, M., P. Axelrad, D. Masters, IEEE Geoscience and Remote Sensing 
Systems (IGARSS) 2001, Sydney, Australia, 9-13 July, 2001. 
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ABSTRACT 

Global Positioning System (GPS) signals reflected from 
the ocean surface have potential use for various remote 
sensing purposes. Some possibilities are measurements of 
surface roughness characteristics from which wave height, 
wind speed and direction could be determined. In this 
paper, recent reflected GPS measurements collected via 
aircraft with a delay mapping GPS receiver, are used to 
explore the possibility of determining ocean surface wind 
speed and direction during Hurricanes Michael and Keith 
in October, 2000. To interpret the GPS data, a theoretical 
model is used to describe the correlation power of the 
reflected GPS signals for different time delays as a 
function of geometrical and environmental parameters. 
Wind direction estimates are based on a multiple satellite 
non-linear least squares solution. The estimated wind 
speed using surface-reflected GPS data collected at a 
variety of wind speed conditions shows an overall 
agreement better than 2 m/s with data obtained from 
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nearby buoy data and independent wind speed 
measurements derived from TOPEX/Poseidon, ERS, and 
QuikSCAT observations. Wind direction agreement with 
QuikSCAT measurements appear to be at the 30 degree 
level. In addition to the actual results we also present an 
error analysis describing the expected performance of the 
GPS technique. 

INTRODUCTION 


The use of GPS as a forwardscatter remote sensing tool 
has come to fruition in the last few years (Katzberg 
[1996], Garrison et al. [1997, 1998, 1999], Komjathy et 
al. [1998, 1999], Lin et al. [1999], Armatys [2000]). 
NASA researchers Stephen J. Katzberg of LaRC and 
James L. Garrison now at Purdue University have 
determined the properties of the ocean-reflected signal 
and have developed a specialized GPS receiver called the 
Delay-Mapping Receiver (DMR) to measure the reflected 
signals [Garrison et al., 1997]. 

Other investigations of ocean reflected GPS signals are 
being conducted by the Jet Propulsion Laboratory (JPL) 
and the European Space Agency (ESA), focused primarily 
on the application of reflected GPS signal tracking to 
altimetry [Martin-Neira, 1993]. These groups have 
conducted a number of experiments from static locations 
and aircraft, and investigated signals received from a 
spacebome antenna [Lowe et al., 2000a, 2000b]. 

Using the DMR and models to predict the interaction of 
the LI GPS signal at 1575.42 MHz, researchers at NASA 
Langley Research Center, Purdue University and the 
University of Colorado at Boulder have been able to 
estimate wind speeds on the ocean surface with an 
accuracy of about 2 m/s. Results to date have advanced 
the understanding of reflected GPS signals and provided 
direct experimental evidence of their application to ocean 
remote sensing and mapping. 

The model employed in this paper’s estimation process 
was developed by Zavorotny and Voronovich (Z-V) and is 
documented extensively in Clifford et al. [1998], 
Zavorotny and Voronovich [2000] and Komjathy et al. 
[2000]. The Z-V model employs a forwardscatter radar 
equation with the geometric optics limit of the Kirchhoff 
Approximation. For backscattering theory and 
applications, the interested reader is referred to Ulaby et 
al. [1986], The wave spectra used by the Z-V model is 
that from Elfouhaily et al. [1 997]. 


The Z-V model takes the form: 
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where / c j| 2 j is the reflected power for any delay 

bin t and Doppler offset f c ; T i is the integration time in 
seconds; SR is the complex reflectivity of the ocean at LI; 
D is the antenna gain; A is the correlation function of 
the GPS C/A code; S is the Doppler sync function; P is 
the probability density function (PDF) of the surface 
slopes; q is the magnitude of the scattering vector q ; R 0 
is the distance from some point on the surface point to the 
GPS satellite; R is the distance from the GPS receiver to 
some point on the surface, c is the speed of light; f D is the 
Doppler shift at the specular point; f c is the compensation 
frequency or the Doppler offset to some point p ; and 
p is a vector from the specular point to some other point 

on the surface. For our aircraft experiments discussed in 
this paper, equation (1) can be simplified and S set to 
unity. 

INSTRUMENT AND DATA 

The use of GPS in a bistatic radar configuration to 
measure surface properties relies upon our ability to 
extract information from the reflected signal. For standard 
GPS navigation applications, the receiver’s main functions 
are to measure the signal delay from the satellite (the 
pseudorange measurement) and the rate of change of the 
range (the Doppler measurement) (see e.g , Parkinson et 
al. [1996]). Conversely, in our remote sensing application, 
the primary measurement is the received power from a 
reflected signal for a variety of delays and Doppler values. 
The basis of this measurement and its sensitivity to the 
surface conditions is discussed in the following 
paragraphs. 

The Delay Mapping Receiver (DMR) is a software 
configurable GEC Plessey (now Mitel Semiconductor) 
GPS Builder-2 receiver modified to observe reflected 
LHCP signals from two GPS satellites and to record 
correlator power at 10 consecutive half-chip intervals. The 
half-chip intervals are analogous to range bins in a radar 
receiver and are used to isolate power reflecting from a 
specific region on the ocean surface. Signals reflected 
from the ocean surface originate from a glistening zone 
(see Figure la) surrounding a nominal specular reflection 
point. The size and shape of the glistening zone are 
functions of the roughness of the ocean surface. To 
measure the reflected power from this glistening zone, the 
receiver-generated pseudorandom noise codes are delayed 
in time with respect to directly received, line-of-sight 
signals. This isolates power originating from the region 
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surrounding the specular reflection point. The shape of the 
resulting waveform of power-versus-delay is dependent 
upon the roughness of the ocean surface (see Figure lb). 
This roughness is in turn a function of the surface wind 
speed and direction, and therefore provides a means to 
retrieve these geophysical parameters from the GPS 
reflected signal power measurement. 



Ideal Autocorrelation 
Function 



Delay Relative to Specular Reflection Point 
(Code Chips) 

Figure lb. Correlation function shapes for ideal direct 
GPS signal and for reflected signals from rough surfaces. 

METHOD FOR WIND SPEED AND WIND 
DIRECTION RETRIEVAL 

Preprocessing the Data . Before the estimator can make 
use of the reflected GPS data, they must be preprocessed. 
Preprocessing takes place in several steps. First, the noise 
floor is computed for each data set. This is done by 
computing the mean of all the points before the First 
correlation peak of the reflected signal. After computing 
the noise floor, it is subtracted from ail data points. These 
reflected data points are then normalized by dividing by 
the total reflected power. Normalization is necessary to 
remove effects of uncalibrated receiver gains. The total 
reflected power is computed by summing all of the 


correlation measurements for one waveform, essentially 
integrating the correlation waveform. Total reflected 
power is chosen for normalization because it should be 
nearly constant due to conservation of energy. Finally, the 
data is broken into one-minute segments for use by the 
estimator. 

An estimate of the path delay, r , is computed using post- 
processed positions of the satellite and receiver. The 
satellite positions are interpolated from International GPS 
Service’s (IGS) 15-minute precise positions, and the 
receiver positions are interpolated using the receiver’s 
navigation solution. Using these positions and the WGS- 
84 ellipsoid, an estimate of the specular reflecting point 
coordinates on the earth surface is computed. The path 
delay is then estimated from these three positions. 
Because the delay variable computed from the receiver 
and satellite geometry with respect to the WGS-84 
ellipsoid may contain errors, a shifting parameter is 
introduced. A scaling parameter is also introduced which 
compensates for errors in the assumptions made during 
normalization of the measured power. Because the total 
power measured over the 10 delay bins fails to include 
power over the same range of delay as the modeled 
waveform, inclusion of a scale factor accounts for this 
discrepancy during normalization. During preprocessing, 
we also quality-check the data and eliminate outliers by 
computing the mean and standard deviation of the 
reflected power in each delay bin. 

Main Processor. The state for the estimation process 
contains wind speed, wind direction and as an option, path 
delay error estimates, and scaling parameters can also be 
simultaneously estimated. For routine data processing, the 
software is able to estimate path delay errors by aligning 
the waveform leading edges. 

The basis for wind direction determination is that the PDF 
of surface slopes is wider in the direction of the wind. 
This produces an asymmetry in the glistening zone. With 
delay measurements from a single satellite, it is not 
possible to unambiguously identify this asymmetry 
direction because the integration over a delay bin or 
annulus tends to obscure the uneven distribution, creating 
an ambiguity with respect to the asymmetry direction. 
Recovery is possible with multiple satellites under the 
condition that the glistening zones are due to the same 
surface wind conditions. Because the annuli for the two or 
more satellites are not mutually concentric, these 
measurements provide the necessary conditions for 
observing the PDF asymmetry. 

In the latest version of our algorithm, we implemented the 
option to process any number of satellites in a single batch 
to fit the measured to the modeled waveforms using a non- 
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linear least squares algorithm in MATLAB. In the 
algorithm, residuals are minimized using a Nelder-Mean 
direct search method to adjust the state. 

By processing two or three satellites simultaneously, both 
wind speed and direction can be solved. To make the 
multiple satellite estimator code run faster, we created an 
extensive waveform database using combinations of 
receiver height, elevation angle, wind speeds and wind 
directions. 

FLIGHTS 

The NOAA Hurricane Research Division, with the 
assistance of NOAA/Environmental Technology 
Laboratory, installed a DMR GPS receiver in one of the 
WP-3D Hurricane Hunters N42RF u Kermit”. Members of 
NOAA Airborne Operations Center at MacDill AFB 
installed the receiver in August 2000. The first data was 
collected from the outer edges of Hurricane Keith on 
October 1, 2000. The emergence of Hurricane Michael 
presented the first opportunity to traverse the core of a 
tropical cyclone on October 18, 2000. Hurricane Michael 
formed in the Western Atlantic the evening of October 16. 
It reached tropical storm strength the next morning and 
was classified a hurricane that same afternoon. Michael 
increased forward speed the night of October 18. It sped 
northward, made landfall in Newfoundland, Canada on 
October 19, and quickly began losing tropical 
characteristics. At the time of the aircraft penetration, the 
storm was moving at approximately 18 m/s. Figure 2 is a 
map derived from the GPS position information, which 
shows the aircraft flight path from MacDill, Florida into 
the center of Hurricane Michael and subsequent flying in 
and around storm. The cross represents the center of the 
storm (from Hurricane Research Division data) during the 
time the GPS surface reflection data was taken. 
Preliminary results and findings of the data processing 
were presented in Katzberg [2001]. 

For comparison purposes, we used wind speeds from 
TOPEX/Poseidon dual-frequency altimeter, ERS 
altimeter, buoy measurements, flight level wind speed and 
on-board step-frequency microwave radiometer (SFMR) 
data (see e.g., Ulaby et al. [1986]). To compare GPS- 
derived wind direction with an independent ground-truth 
we used QuikSCAT-derived wind measurement. For 
detailed discussion of different satellite wind speed 
retrieval techniques, see Stewart [1985]. 

RESULTS 

Wind Speed Retrieval. The aircraft flew out from the coast 
at an altitude of 4500 meters. It descended to 1400 meters 
traversed the eye of the storm and then descended further 



Figure 2. Flight trajectory for Hurricane Michael of 
October 18, 2000 adapted from Katzberg [2001]. 


to 500 meters where it stayed for the bulk of the time the 
GPS equipment was operated. At approximately 15:50 
UT, the aircraft flight path crossed a TOPEX ground 
track. Figure 3 shows the GPS-derived wind estimates at 
this time, for each satellite separately. A combined 
solution using PRNs 15, 21 and 23 is presented in Figure 
4. In Figure 3, we also plotted the mean of the individual 
solutions along with the standard deviations. The 
estimated wind speed estimates using single satellites 
ranges between 6 m/s and 10 m/s. The combined solution 
in Figure 4, agrees with the TOPEX solution within an 
RMS of 0.7 m/s. The RMS of the individual satellite 
estimates is 1 .2 m/s. 

At approximately 14:55 UT, the aircraft crossed the ERS 
ground track. We again processed the satellites separately 
as shown in Figure 5 along with the mean and standard 
deviation. A combined solution was obtained showing 
better agreement with the corresponding ERS ground- 
truth (0.7 m/s RMS, see Figure 6) than the average 
solution from individual satellites (0.9 m/s RMS). We 
believe that the larger differences near the start of the 
graph arc the result of a larger separation between the 
aircraft position and the ERS ground track at that point. 



lime in HI 

* PRN2I « PRN 1 5 * PRN23 PRN29 * PRN3Q -•-TOPEX -*“Mean 


Figure 3. Satellite by satellite GPS wind speed (WS) 
solution for TOPEX pass. 
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time In LIT 

PRN2 1- 1 5-23 — TOPHX 


Figure 4. Multiple satellite GPS wind speed (WS) solution 
for TOPEX pass. 



• PRN2I ■ PRN29 * PRN09 PRMI5 * PRN17 • PRN23 -m- ERS —•—Mean 

Figure 5. Satellite by satellite GPS wind speed (WS) 
solution for ERS pass. 



14:52 14:53 14:54 14:55 1456 1457 14 5» 1459 I5KK) 15:01 
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Figure 6. Muliple satellite GPS wind speed (WS) solution 
for ERS pass. 

In Figure 7, we show the combined solution for a three- 
hour segment of the flight track that also included the 
TOPEX and ERS data segments. Also displayed is the 
aircraft altitude. Time tags corresponding to TOPEX and 
ERS passes, and buoy measurements are superimposed 
showing generally a good agreement between GPS wind 
estimates and ground-truth for the pre-hurricane part of 
the flight. Furthermore, data from Hurricane Research 
Division of NOAA representing “flight level wind speed” 
(FLWS) was obtained and plotted in Figure 7. FLWS data 
was provided in 1 -minute and 10-second time series. In 
Figure 7, due to similarities between the two time series, 


we only plotted the 1 -minute averaged data set. FLWS 
data is derived from the drift the aircraft experiences when 
flying into the storm. Also displayed is the time series 
from step-frequency microwave radiometer (SFMR). At 
approximately 16:30 UT, the aircraft was approaching the 
hurricane region. This is apparent in rapid increases in the 
wind speed from both FLWS and SFMR data. At around 
16:55 UT, the aircraft flew over the eye of the hurricane. 
The high variability of the GPS-retrieved wind speed 
made the comparison difficult. The peak GPS wind speeds 
seem to be consistent with the SFMR data, but capturing 
the eye of the hurricane was unsuccessful. It is possible 
that this is because the rapid movement of the storm did 
not allow for the condition of fully developed seas 
assumed in the wave spectrum and Z-V models. 
Furthermore, the aircraft descending from 4500 to 1400 
meter altitude made waveforms become sharp and narrow 
resulting in the modeling and retrieval process more likely 
to fail. 

In Figure 8, we show an example for 15:44 UT (see 
Figures 3 and 4) of the estimator convergence in terms of 
both the wind speed estimate and the sum of the squared 
measurement residuals. The filter was initialized with a 
wind speed of 10 m/s and scaling parameter set to unity. It 
is shown that after approximately 70 iterations, the wind 
speed converges to 7.8 m/s. In most cases, convergence is 
reached in fewer than 30 iterations when the state 
estimates are initialized with the solution from the 
previous segment wind speed estimate. 

Wind Direction Retrieval. The second data set we present 
is for October 1, 2000, taken from the outer edges of 
Hurricane Keith in the Gulf of Mexico to the west of 
Florida. Wind speeds ranged from 6 m/s to 10 m/s. What 
distinguishes this data set from Michael is the fact that we 
have QuikSCAT wind field data available that was taken 
within one hour of the GPS measurements. In this case, 
we retrieved wind speed as well as wind direction 
information using the multiple satellite solution from 
PRNs 01, 03 and 13 as described earlier. 

In Figure 9, we plot the GPS-derived wind speed and 
direction solution superimposed on the QuikSCAT wind 
field plot. The aircraft was flying from west to east 
direction. Wind speeds are plotted next to the base of the 
arrows representing the GPS-estimated wind directions. 
The left part of the flight track shows a good agreement in 
both wind speeds and wind directions. 
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Time in UT 

GPS Estimates FLWS(lmm) FLWS (10 sec) — * — SFMR 

Figure 7. GPS wind speed (WS) estimates along the flight path for Hurricane Michael of October 1 8, 2000. 




Figure 8. An example of wind speed (WS) convergence 
showing dependence between number of iterations, 
estimated wind speed and residuals for the TOPEX pass. 



Figure 9. GPS-derived wind vector estimates for the 
vicinity of Hurricane Keith of October 1, overlaid on 
QuikSCAT wind field measurements. 


The flight path on the right hand side indicates fluctuating 
wind direction with stable wind speeds. We believe that 
the reason for the unstable wind direction estimates is 
rapidly changing aircraft altitude as it was getting closer to 
the Florida coast. The first few GPS wind direction 
estimates also show some fluctuations even though the 
wind speeds are consistent with QuikSCAT data. The 
discrepancy in wind direction might be due to the fact that 
our estimator flagged too many data points as outliers 
resulting in smaller number of waveforms to be available 
to the estimator. 


It is interesting to see that GPS-derived wind is mostly in 
south-west direction, whereas QuikSCAT data shows 
similar wind direction only above 24°30\ Below this 
latitude QuikSCAT wind is predominantly westerly. One 
possible explanation is that our bistatic GPS scattering 
model assumes well developed seas. This simplifying 
assumption could limit our accuracy in this coastal region 
where the Cuban land mass is approximately 50 km to the 
south. Furthermore, ocean currents typical for the Cuban 
coastal region, not taken into account in our model, may 
also have an impact on our GPS-derived wind field. The 
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fact that we have good agreement with QuikSCAT in well 
developed sea conditions points us to the need for further 
studying the effect of un-developed seas conditions. 


In Figure 10, we plotted an example for estimated wind 
speed, wind direction and residuals to demonstrate the 
convergence of the solution on the accepted truth. It is 
shown that after 40 iterations the solution converges on 
7.6 m/s wind speed and about 30 degree wind direction. 
QuikSCAT indicates 7.6 to 7.8 m/s wind speed and about 
40 degrees wind direction. 
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Figure 10. An example for wind speed (WS) and wind 
direction (WD) convergence for October 1 , 2000. 


ERROR ANALYSIS 


the residuals occurs for wind speeds between 6 and 8 m/s. 
On the other hand, in Figure 12, we show that the 
minimum of the residuals are obtained for wind direction 
between 20 and 40 degrees. The individual curves 
represent a solution with different wind speeds between 3 
and 15 m/s starting with 3 m/s at the top. We can clearly 
see that the difference between measured and modeled 
values (the residuals) are at the minimum when the right 
answer (WS and WD) is used to compute the 
measurement residuals. In Figure 12, we can realize that 
in the vicinity of the QuikSCAT result, (7 or 8 m/s) there 
exist solutions with combinations of wind speeds and 
wind directions resulting in the same sum of the squared 
residuals. To demonstrate this, we plotted the residuals 
using the combinations of wind speeds and wind 
directions. In Figure 13, iso-residual lines show possible 
combinations of wind speeds and wind directions resulting 
in the same minimum of the sum of the squared residuals. 
The “hole” in the middle shows a set of solutions close to 
the final solution. However, the estimator indicates that 
the sum of the squared residuals are at the minimum at the 
middle of the “hole” that can be characterized with wind 
speed of 7.6 m/s and wind direction of 30 degrees which 
is in very good agreement with the QuikSCAT 
measurements. 

A reliable wind direction estimate cannot be obtained 
using delay measurements from a single satellite. This is 
demonstrated in Figure 14. Although the estimated wind 
speed is about 8 m/s, the wind direction is estimated to be 
near 140 degrees, coinciding with the incidence angle of 
the satellite with an additional ambiguity at about 20 
degrees. There are two types of ambiguities in the GPS 
wind direction retrieval. Firstly, there is a 180 degree 
ambiguity, related to the symmetry of the PDF of slope in 
up/down wind direction. This ambiguity cannot be 
resolved even using multiple satellites. Secondly, there is 
the ambiguity related to the symmetry with respect to the 
incidence plane. Indeed, there are two possible angles 
between the incidence plane and wind direction that create 
the same GPS reflected signal power. Therefore, the 
ambiguity can only be solved for using additional 
simultaneously observed satellites. In Figure 15, we 
plotted the map of wind speeds and wind directions using 
one satellite only. The estimator shows a larger pool of 
possible solutions before finally arriving at 8 m/s wind 
speed and 170 degree wind direction (aligned with the 
incident plane). 


In this section, we investigate the non-linear wind speed 
and wind direction estimator convergence properties. In 
Figures 11 and 12, we computed the combined solution 
residuals using all possible combinations of wind speed (1 
m/s increments) and wind directions (10 degree 
increments). In Figure 1 1, it is shown that the minimum of 


We also computed the repeatibi lity of the solutions based 
upon consecutive and independent 1 -minute data 
segments with the same three satellites. We assume that 
the wind speed and wind direction do not change over a 
10-minute window. Computing the repeatibility gives a 
measure of the effect of measurement noise on the actual 


Presented at the ION 2001 National Technical Meeting, Long Beach, CA, January 22-24 


1 


solution. This analysis shows that over a 10-minute arc, 
the standard deviations of the wind speed and direction 
are 0.7 m/s and 9 degrees, respectively. As to the accuracy 
of our GPS-derived wind direction estimates, the 
processed data using the combined solution indicated a 
better than 30 degree agreement with the QuikSCAT 
measurement. This is encouraging since the overall 
reported QuikSCAT wind direction accuracy is about 20 
degrees. 



WS rn m's 


Figure 11. Residuals (unitless) versus wind speed (WS) 
for October 1 data set. 



WD in Degrees 


Figure 12. Residuals (unitless) versus wind direction 
(WD) for October 1 data set. 



WD in Degrees 


Figure 13. Wind direction (WD), wind speed (WS) and 
residuals map for the combined multi-satellite solution. 



WD in Degrees 


Figure 14. Illustration of unreliable WD estimation using 
a single satellite (PRN03). 



WD in Degrees 

Figure 15. Wind direction (WD), wind speed (WS) and 
residuals map using a single satellite PRN03. 
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CONCLUSIONS AND FUTURE RESEARCH 

We have demonstrated wind speed and wind direction 
retrievals using a novel multi-satellite approach combined 
with non-linear least squares estimation. GPS-derived 
wind speed and wind direction is compared wind TOPEX, 
ERS, buoy and QuikSCAT measurements. We found that 
processing surface reflected GPS signals in a combined 
solution, as opposed to a satellite by satellite solution, 
gives us better wind speed agreement, at the level of 2 
m/s, with other independent techniques for a variety of 
wind speed conditions ranging from 3 to 10 m/s. A 
comparison between GPS-derived wind direction and 
QuikSCAT wind field showed a better than 30 degree 
agreement in wind direction. We also demonstrated that it 
is not possible to estimate wind direction with delay 
measurements from a single satellite. 

Future research will investigate the performance of the 
technique during high (> 1 5 m/s) wind conditions and in 
coastal regions. Additional model enhancements to 
account for more variable conditions and improvements to 
data quality and processing are required to achieve 
progress in these more challenging conditions. 
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ABSTRACT - Previous papers have demonstrated the 
ability to determine wind speed and direction from ocean- 
reflected GPS signals observed at aircraft altitudes to 
within 2 m/s and 20 degrees of in situ measurements. 
This paper presents a preliminary step towards extending 
GPS-based wind vector retrievals to spacecraft altitudes. 
The work takes the form of visibility, sensitivity, and noise 
analyses. The visibility analysis reveals the number of 
ocean-reflected GPS signals visible per day, a value that 
drives the uncertainty in wind vector estimates. The sen- 
sitivity analysis determines which delay /Doppler cells of 
a glistening zone visible from space are most sensitive to 
changes in wind speed and direction. Finally, the noise 
analysis uses results from the previous two analyses to 
determine the uncertainty in wind vector retrievals using 
the cells identified in the sensitivity analysis. 

INTRODUCTION 

The GPS has been used as a remote sensing tool since 
1997 [1]. Researchers at the University of Colorado (CU) 
and NASA-Langley Research Center (LaRC) developed 
algorithms for estimating wind speed from ocean-reflected 
GPS signals collected at aircraft altitudes. In their work, 
the CU and LaRC researchers demonstrated wind speed 
estimation to within 2 m/s of available in situ measure- 
ments [2]. In 1999, the authors developed algorithms for 
estimating both wind speed and direction from reflected 
GPS signals [3]. The wind speed and direction estimates 
agree with QuikSCAT measurements to within 1 m/s and 
10° under ideal conditions. Komjathy further developed 
the algorithms [4]. 

Spaceborne measurements of ocean-reflected GPS signals 
are the next obvious step in development of GPS-based 
remote sensing. In the summer of 2000, CU teamed with 
the University of Texas at Austin, Purdue University, the 
Applied Physics Laboratory at Johns Hopkins Univer- 
sity, the Jet Propulsion Laboratory, the National Oceanic 
and Atmospheric Administration, and Ball Aerospace to 
propose a student-built satellite to make both wind vec- 
tor retrieval and altimetry measurements [5]. The LEO 
satellite, named SuRGE for Student Reflected GPS Ex- 
periment, would be the first satellite dedicated solely to 
measuring reflected GPS signals. This paper describes 
preliminary work to extend wind retrieval algori thins to 
spacecraft altitudes. 


* Submitted to IGARSS 2001, Sydney, Australia, 9-13 July, 
2001 


VISIBILITY ANALYSIS 

Unlike a conventional scatterometer , the locations of mea- 
surements made by a bistatic radar are not controlled by 
the instrument, but rather by the geometry of the trans- 
mitter, receiver, and surface. To evaluate the Earth cov- 
erage for GPS observations from SuRGE we performed 
a visibility analysis using Satellite Toolkit® (STK) from 
Analytical Graphics and our own MATLAB® scripts. 
STK® was used to produce position and velocities of the 
GPS satellites and SuRGE at 1 s intervals for a 1 day pe- 
riod. A MATLAB® script then computed the positions 
of specular points for all satellites using the iterative algo- 
rithm proposed by Wu, Meehan, and Young [6]. Antenna 
patterns for SuRGE were applied to mask out specular 
points outside the field of view of the reflection antennas, 
and a low-resolution land mass model from MATLAB® 
was also applied to mask out specular points over land. 
The remaining list of specular points represent the avail- 
able bistatic ocean observations over one day. Raw ob- 
servations in the receiver are produced at an integration 
time on the order of 0.3 ms. (This time is selected to corre- 
spond to the coherence time of the area of interest around 
the specular point). Raw observations are averaged over 
1 s or longer to reduce the effects of noise. To provide a 
comparison to QuikSCAT, we have grouped the observa- 
tions based on a surface grid. The grid divides the Earth 
surface into cells of equal area (625 km 2 ). The number 
of consecutive raw observations in each cell is computed. 
This determines the number of samples that can be av- 
eraged to form a wind estimate for the cell. A single 
wind measurement includes all the consecutive observa- 
tions within a specific cell. Thus all measurements cover 
the same surface area on the earth. More observations 
within the cell produce more accurate measurements as 
will be described in the following sections. 

Fig. 1 shows a histogram of the visibility results cate- 
gorized by the number of observations per measurement. 
The bars represent measurements including 3000, 7000, 
10000, etc 0.3 ms observations. Equivalently this corre- 
sponds to averaging times of 0.9 s, 2.1 s, and 3 s. For 
example, the third bar on the histogram shows that al- 
most 7000 wind measurements on a 625 km 2 area can be 
made with an averaging time of 2.1 s. A total of 35,000 
measurements can be made with a minimum averaging 
time of 1 s. 

MODEL EQUATION 

The model employed in the noise analysis was developed 
by Zavorotny and Voronovich (Z-V) and is documented 


extensively in [7]. The Z-V model employs a forward- 
scatter radar equation with the geometric optics limit of 
the Kirchhoff Approximation. The wave spectra used by 
the Z-V model is that from Elfouhaily, Chapron, Katsaros, 
and Vandemark [8]. 

SENSITIVITY ANALYSIS 

The sensitivity analysis revealed which delay /Doppler 
cells on the Earth’s surface show the most sensitivity to 
changes in wind speed and direction. Model waveforms 
were generated at 1/2 C/ A code chip and 500 Hz Doppler 
increments and a gradient analysis was performed be- 
tween adjacent waveforms to determine those cells with 
the largest gradient. The gradient analysis reveals that 
delay /Doppler cells near the specular point provide the 
most sensitivity to changes in wind speed. These cells 
should be used for wind speed retrievals from space. Fur- 
thermore, for low wind speeds the delay /Doppler cells fur- 
thest from the specular point also provide some sensitivity 
to changing wind speed and direction. 

Fig. 2 shows a plot of percent change in reflected signal 
power for a 1 m/s change in wind speed at three represen- 
tative wind speeds: 4, 10, and 20 m/s. Power is plotted 
only as a function of delay in C/A code chips rather than 
versus delay and Doppler. Fig. 2 demonstrates that the 
cells nearest the specular point corresponding to delays 
near zero provide the most sensitivity to changes in wind 
speed. Furthermore, higher relative power values in the 
4 m/s curve show that changes in wind speed are more 
easily detected at lower wind speeds. Fig. 2 also shows 
the two regions of sensitivity for glistening zones caused 
by lower wind speeds. The reason for the two maxima in 
the 4 m/s curve is visible in Fig. 3. Lower wind speeds 
produce cross-correlation functions shown in Fig. 3 with 
larger slopes than for higher wind speeds. The rate of 
change in slope in the tail of the auto-correlation function 
causes the second maxima at low wind speeds. 

NOISE ANALYSIS 

The noise analysis determines the uncertainty in wind 
power measurements from space, measurements which 
drive the wind speed and direction estimates. The noise 
analysis focuses on those delay/Doppler cells identified in 
the previous section and uses typical numbers of observa- 
tions. The following equation from Ulaby describes the 
uncertainty in the power measurements [9]: 


’ N 


( 1 + ^r) + (sr) 

N r 
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adjacent modeled waveforms. The sensitivity analysis in 
the previous section showed that the delay /Doppler cells 
nearest the specular point and near the edge of the glisten- 
ing zone showed the most sensitivity to changes in wind 
speed Using the waveforms corresponding to 4, 10, and 
20 m/s, uncertainties for the all delay /Doppler cells were 
computed at each of the number of observations repre- 
sented on the x-axis in Fig. 1. 

Fig. 4 shows the uncertainty in wind speed estimates 
versus C/A code delay as a function of number of obser- 
vations. The lines in Fig. 4 correspond to the different 
numbers of observations presented in Fig. 1 and are con- 
stant because they were computed as the ratio of uncer- 
tainty or standard deviation in power per delay bin to the 
power for each delay bin. From (1) it is evident that this 
ratio is a constant for all delay bins. The set of curves 
in Fig. 4 is identical for all wind speeds because of the 
constant ratio as well. 

Using (1) and Fig. 2, we can estimate the wind accuracy 
as a function of number of measurements in the estimate. 
Dividing (1), we find the ratio of the power measurement 
error to the total power. For 3000 observations, the power 
measurement error is 4%. Based on Fig. 2, we can see 
that for a wind speed of 4 m/s a power measurement near 
the specular point with an uncertainty of 4% corresponds 
to a wind uncertainty of approximately 0-25 m/s. For 
10 m/s wind the same power measurement error gives a 
1 m/s wind uncertainty, and about 4 m/s for a 20 m/s 
wind. 


CONCLUSIONS 

This paper presents several key preliminary results re- 
lated to estimating wind speed and direction from space 
using reflected GPS signals. First, predicted numbers of 
reflected GPS signals from space for 625 km 2 equal-area 
cells are presented. Second, delay/Doppler cells near the 
specular point are identified as having the greatest sensi- 
tivity to changes in wind speed. It is proposed to perform 
a similar analysis for identifying wind direction sensitiv- 
ity. Finally, uncertainties in reflected signal power predict 
wind speed retrievals to within 1 m/s at low wind speeds 
and to within 2 m/s at high wind speeds for all equal-area 
cells providing at least 1 s of averaging time. 
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where a 9 is the standard deviation of the signal, P s is 
average power of the signal, S n is signal to noise ratio 
for a single 0.3 ms observation, and N r is the number of 
observations. In our simulations, we assume that S n = 1, 
a pessimistic scenario. 

Uncertainties in wind speed and direction measurements 
can be derived using the above equation and comparing 
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Figure 2: Plot of the percent change in reflected signal 
power for 1 m/s change in wind speed versus delay. Curves 
correspond to 4, 10, and 20 m/s. 



Figure 3: Plot of power versus delay for wind speeds of 3, 
5, 10, 14, and 20 m/s. 
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Figure 1: Histogram of number of measurements made in 
1 day formed from 3000, 7000, 10000, etc. raw observa- 
tions within a 625 km 2 surface cell. 
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Figure 4: Plot of the percent uncertainty in reflected sig- 
nal power as a function of delay for 4 m/s and 3000 obser- 
vations. Wind speeds of 10 and 20 m/s produce identical 
plots. 
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Overview 


Multipath is a common error source affecting GPS use for positioning and navigation. However, 
Delay Mapping researchers have only recently recognized, that a special kind of multipath - that from GPS signals 
Receiver reflected off the sea surface - can be used as a tool for oceanographic remote sensing. 



Because ocean surface roughness affects GPS signal reflection, we use the resulting multipath signal 
to determine such factors as wave height, wind speed, and wind direction. The reflected signal's 
strength is also a discriminator between wet and dry ground areas and, therefore, could be applied to 
coastal and wetland mapping. We are also looking at the reflected signature from sea ice. 


The Colorado Center for Astrodynamics Research at the University of Colorado, collaborating with 
researchers at NASA, is currently experimenting with data from software configurable receivers flown on 
aircraft platforms to retrieve ocean surface wind speed. In conjunction, modeling efforts by scientists at 
the NOAA/Environmental Technology Laboratory in Boulder support this research. 

An introduction to the topic was published in the April 1999, issue of GPS World . A concise overview is 
also presented by Dr. Valery Zavorotny. 
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